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We investigate a system of dense polyelectrolytes in solution. The Langevin dynamics of the 
system with linearized hydrodynamics is formulated in the functional integral formalism and a 
transformation made to collective coordinates. Within a dynamical Random Phase Approximation 
(RPA) integration over the counter- and salt ions produces the Debye-Hiickel-like screening of the 
Coulomb interactions with dependence on the frequency only as part of a more complicated coupling 
structure. We investigate the dynamics of the structure factor as well as the collective diffusion 
coefficient and comment upon the viscosity of the whole system of polymers with counterions and 
fluid in the simplest approximation. The coupling of the various components of the system produces 
nontrivial diffusive behavior. We draw conclusions about the relationship of the three length scales 
in the present system, i.e. the static screening length, the hydrodynamic screening length and the 
Debye length. 

I. INTRODUCTION 

The statistical mechanics of polyelectrolytes has recently been the object of much research, of both analytical 
nature and simulations, ranging from a description of the flexible single-chain propertieslttil to the properties of the 
dense melts in collective coordinates (refs 8-10 and references therein). The presence of the long-ranged Coulomb 
interactions makes such systems completely different to uncharged ones. Underlying many of these approaches is the 
electrostatic screened interaction according to Debye-Hiickel, of the form 

VDu(r) = k B T^e- K ^. (1.1) 

Here k represents the screening parameter, or its inverse the Debye-Hiickel screening length, and As is called the 
Bjerrum length, 

Here e is the dielectric constant of the fluid, and fcgT the temperature multiplied by the Boltzmann constant. The 
parameter, k, is proportional to the Bjerrum length and to the concentration of the counterions. Its inverse is tha 
Debye screening length. The effective potential can be derived from the linearized Poisson-Boltzmann equationtil 
or by utilizing collective density variables together with the random phase approximation in statistical mechanics 
(RPA) approachtl and integration over salt and counterion degrees of freedom, which procedure again repipduces that 
form of effective screened interaction. In other models the mean spherical approximation (MSA) is used.t2l However, 
whether the Yukawa form of the pairwise potential is truly accurate for macromolecular systems still remains a 
debatable question and even for the understanding of the critical behavior in thc-casc of the simple electrolyte a 
debate concerning an appropriate minimal model exists (See, e.g. Levin and FisherEf). 

On the statistical mechanical level, the properties of the single chain have been under much discussion as to the 
precise dependence upon k of the properties of the chain, which is expected to behave as with an enhanced persisted 
length due to the electrostatic interaction. Recent approaches-to the problem are a treatment by Barrat and Joannym 
and a renormalization treatment by Liverpool and |S±apper£l On the level of systems also of higher than extremely 
dilute concentrations recent work by Muthukumanlil has been based upon the assumption of an electrostatically 
altered effective, Gaussian step-length. As regards the implementation of the RPA it has also been shown! that this 
approximation produces qualitatively correct results. It needs to be altered when strong, attractive interactions can 
cause complexationlla or when the effects of the strength of the potential are such that the Gaussian assumption must 

The dynamics of polyelectrolytes has recently been addressed in scaling-theoretic_approach by Dobrynin et a/0'0 
An important experimental result in those dynamical considerations is the Fuoss lawlia which states the proportionality 
of the viscosity and the inverse square root of the n n ly' nTI concentration at sufficiently high densities. This was recently 
confirmed by Dobrynin, Rubinstein and ColbyEj'O P-ttfr treatments of this law have been by Witten and PincusE£l 
and by RabinnJ Previously, a theoretical treatmentOEl of the polyelectrolyte solution viscosity has been given by 
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one of the authors by making use of the results of a mode- mode coupling approximation by Hess and KleinEa inpwhich 
a peak of the viscosity was derived. Recent experimental work by Isambert, Ajdari, Mitnuk, Viovy and Prostcfl has 
highlighted a physical instability of a polyelectrolyte solution occurring when this solution is subjected to an external 
electric field. 

In this paper we formulate the dynamics of a system of polyions with salt and counterions in the Langevin approach, 
and couple the dynamics of these objects to a solvent which is assumed to be governed by linearized hydrodynamics. 
It is possible to bring the different length scales inherent to this problem, i.e. the charge screening length, the 
hydrodynamic screening and the polymer length, into relation with each other. It is assumed that the system is 
sufficiently dense and the chain conformations close to the Gaussian ideal to permit the use of a random phase 
approximation (RPA) in the context of the dynamics (see, e.g. ref 25). Furthermore, it is assumed that linearized 
hydrodynamic properties of the solvent apply and that the investigation is limited to the semi-dilute regime in which 
the effects of entanglements of the macromolecular strands can be neglected. Importantly we find that this coupling 
of the dynamics produces results which are at variance with calculations where an effective Debye-screened potential 
were used straight away. 

The paper is outlined as follows. Section II introduces the dynamical equations of the model system. These 
expressions are xeformulated in terms of a set of -collective variables and expressed through the formalism of Martin, 
Siggia and RoseE£l and of Peliti and De Dominicis.EZl The subsequent sections and subsections under III show the results 
for computations of the correlation functions: this is the dynamical structure factor from which we draw conclusions 
about the diffusional behavior. The final section recounts the results, discusses the hydrodynamic aspects, and 
concludes on the interplay of the screening lengths. 

II. THE MODEL 

A. The Langevin Equations 

The model consists of the Edwards Hamiltonian for the flexible polyions with the excluded volume interaction 
written up to second virial order with the electrostatic potential containing the appropriate Bjerrum length, Xb = 
e 2 /(ATreksT), for the solvent. The dielectric constant for the solvent is given by e = e r eo, k B T the Boltzmann 
constant multiplied by the temperature, and e is the electronic charge (For the sake of simplicity only monovalent 
systems are discussed, although changing valences would present no significant problem in the formalism as used and 
approximated here.). Furthermore, counterions and salt ions occupy the solvent such that overall electrical neutrality 
is maintained. We label the four different species of molecule, i.e. polyion, counterion, and the two types of salt ion, by 
Greek subscripts a = 1, . . . , 4. There are exactly N a identical ions of each species, with N3 = N4, which have position 
coordinates given by r pQ (s, i), where the contour length s is only a relevant variable for the macroions (a = 1). Apart 
from this no further internal degrees of freedom are associated with any of the system constituents. The complete 
Hamiltonian expresses the flexibility of the chains, the excluded volume contributions and the (unscreened) Coulomb 
interaction. 

_ 3^ t L ' f dr pl ( S ,t) \ 2 1^^^ ffi [ L ° ft 

H = Ye^ ds {—d^J + 2"EEE E J I ds Pa ds g ,x 

p=l JU V 7 a = l 13=1 p a = l qi3=l J ,/U 

x ( w aP 5 {v pa {s la , t) - r qP ( Sj , , t)) + M a0 (1 ~ (1 - biA ) (2.1) 

Y \Ypa\Si al t) r g/H S j/3 > V I J 

The Kuhn length of the chains is given by £, and the numerator which contains Kronecker delta functions of the 
term expressing the Coulomb interaction is a technical description to exclude self- interaction by the point particles. 
Formally, in order for the general form above to make sense, the lengths, L ai of the counterions and salt ions are 
one monomer. The Coulomb interaction is represented by the matrix M a p of the appropriate interaction between the 
polymer and the counterions of the system: 

/ P -f +f -f \ 
M aP = X B k B T ~f +J ~\ +H . (2.2) 

W +1 -1 +1/ 

The condition of electrical neutrality implies that N%Lf — N2 + N3 — N4 = must hold. The quantity / represents the 
constant, linear charge density of the polymer. The matrix w a p represents the excluded volume interaction between 
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the various components of the system. The Hamiltonian and the hydrodynamics then lead to the coupled, stochastic 
Langevin equations of the system for the chain molecules, the counterions, and the fluid velocity field.EH When stick 
boundary conditions are employed the stochastic equations can be written as 

= -J^paM +v[r pa (t)] - IJ^+^t) = (2.3) 



4 iV Q Jff 



-EE^7t) 5(x_rpa( * ))+f,;(x ' i) r = a (2 - 4) 

a=2 p=l P a \ ) ) 

The random forces, f pa (spaji) and f„(x, i) for the dynamics of the molecules and the fluid dynamics, respec- 
tively, are correlated as Gaussians and satisfy the fluctuation-dissipation theorem given by (f pa (s Pa , t)f q p{s qi} , i')) = 
-£-5 a pl5 Paq p5(t — t')5(s Pa — s q0 ) and (L,{x,i}f v (x' ,t')) = — 2r]l\7- 2 .5(x — x')S(t — t'). We are also working in the case 
of the incompressible fluid continuumES and, hence, impose the transverse condition upon the velocity field of the 
fluid of bare viscosity ij. The components have friction coefficients £ Q which are in general different. The energy units 
above have been chosen such that ksT = 1 and for subsequent sections the distance units of the polymer will also be 
chosen such that the prefactor 3/2£ in the Hamiltonian vanishes. 



B. Functional Integral Formalism 



Wg-*iow treat the set of coupled differential equations by rising the functional integral approach of Martin, Siggia and 
RoseEa and follow the scheme of Fredrickson and HelfandcS used for uncharged solutions of monodisperse polymers. 
Here we shall be concerned primarily with the collective behavior of the system dynamics, and, in following the work 
of Fredrickson and Helfand, we shall later ignore contributions in the collective variables which are of order higher 
than two. 



By making use of causality and that the Jacobian of the transformation then becomes unitycrl the functional integral 
expressions become: 



Z 



4 N a 



IJ Y[^r pa {t)Vr pa {t) \ Dv(x,i)X*(x,t) 

4 N a „ +00 pL a 



exp < +i 



EE 



dt 



dsr pa {t) ■ Lp 



+ 00 



dt / d xv(x,t) ■ L v 



(2.5) 



Angular brackets indicate the averaging over Gaussian noise and £1 represents the volume. Following ref 25 we 
introduce collective variables in terms of collective density and response functions for each species a (with a number, 
N a of molecules) : 



p a (k,t) = E / ds exp(+ik-r pa (s,t)) 

p=l J( > 

er Q (k,i)=^ / dsr pa (s,t) exp(+ik-r pct (s,t)) 



P =i J o 



JVi 



P =l 



ds 



ds 2 



exp (+zk • r p i(s,t)) . 



(2.6) 
(2.7) 
(2.8) 



Clearly, the 7r a has no meaning when it relates to point-like salt or counterions in the system. 

If one were to impose the additional condition that the whole system finds itself in a spatially varying, external field, 
which oscillates at a frequency Q, such as 2V(k) cos Hit the Langevin equations in the MSR formalism of equation 2.5 
would have an additional contribution in the exponent of: 
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£ ex t = EE 7^^( k ) k ' (^(k, w + w) + o- Q (k, w - w)) . 

a u k 



(2.9) 



The linearized hydro-dynamic formulation used here permits the integration over the fluid velocity fields to be per- 
formed exactly and enables the determination of the frequency-dependent structure factor of the molecular inclusions 
from which the diffusion coefficient of these components is derived. Consequently, one can write 



Z = J {jjDr£>r}exp-£, 



(2.10) 



where £ = Cq + C\ + £3. The variables and integrals are all expressed in terms of their of Fourier representation, the 
random forces have been integrated and the functional integral over the velocity fields, which occur only to quadratic 
order, have been performed. The functions above are given by 



£0 




1 d 2 

iojR pl (s,u>) + — T —rR pl (s,uj) 
Ci c* s 



-Rpi(s,u) 



£1 = 



+EE f p Q ( w ) 

a=2p=l 
4 4 

* E E [ [ ~^ pa (~ k > ^ [ 4?rfc 



7-fpaM 

Set 



(2-11) 



a=l 0=1 " 
4 4 

+EE 

a=l 0=1 

rjk 2 



2 M a p + w a p\ ik ■ erg(k,w) 



Ca = 



LU 2 + rj 2 k4 

444 

*EEE 

a=l 0=1 7=1 
\-2 



(l-kk)-<r a (-k,-w) 

-kk)-cr /5 (k,^) + — 

ijj -\ 

1 



irjk 



r 7T /3 (k,w) 



kk' 



+ irjk 2 



(l-kk)-^(-k,-a;)-k'p/3(k',w') 



x ^^(k') -2 + w Pl ] p 7 (k-k',w-w') 



(2.12) 



(2.13) 



The friction coefficients for the polyion chain segments and the counterions are given by d and £2, respectively. 
The second term of £7 and the contribution £3 come from the coupling to the hydrodynamic medium with sticking 
boundary conditions. 

Since we assume to be dealing with sufficiently dense systems below the onset of entanglement effects, we can make 
use of the approximation that the chains behave as Gaussians. The random phase approximation is implemented 
utilizing exp— Cq as the statistical measure for the collective variables as done by Fredrickson and HelfandEa The 
results are summarized in the appendix. Importantly, in the quadratic approximation the collective variables p a and 
its dynamical conjugate field k ■ er a decouple from (1 — kk) ■ er Q and (1 — kk) ■ n a . This decoupling and the quadratic 
distribution also cause the cubic contribution £3 to be a negligible contribution such that the hydrodynamics will not 
have any effect upon p a and k • cr a when the interactions are considered at this order. 

For the sake of simplicity, we perform the calculation for the system consisting of polyions and the counterions 
leaving the component labels a — 1 and a — 2 with concentrations given by pi and P2, respectively. Extending the 
results for salt ions is straightforward. Investigating the properties of the structure factor requires the use of the 
collective variables which can be contained in a supervector, 



^(k,w) = (pi(k,w), k-£ri(k,w), P2(k,v), k-<7 2 (k,cj)). 
Averages concerning these quantities can now be expressed as follows: 



(A) 



RPA 



N Vi/jAexp- V T (-k,w) -T- i>{k,uj) 



(2.14) 



(2.15) 



Where the matrix T contains the contributions from the RPA and the interaction term, £/. An explicit form of this 
matrix is given in the following section. 
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III. THE DYNAMIC STRUCTURE FACTOR 



After including the interactions and the two coupled density fields which are not in the transverse direction, the 
computation of the structure factor reduces to the identification of the appropriate elements of the inverse of the 
matrix T given below. The excluded volume interaction is assumed only between polyion chain segments, such that 
we take the matrix w a p to be given by the expression w a p = 5\ a 5i(3V. Here we consider the case of large length scales 
and large times. These conditions, kR g <C 1 and Lur r <C 1, where in the units which have been chosen R g = L/2, and 
the Rouse time, r r 



L Ci7r , lead to the matrix in the quadratic supervector formulation, 



i 

Lpi 



pik 2 



i (lLO 
h — f- 

Lp 1 p x k 2 
iXsf Ci 



k 2 



k 2 



V 





i*B f 
k 2 



k 2 pi 
i^B f 
k 2 





% 

P2 





k 2 


(20J 



1 

P2 



k 2 


C2W 



Pi 



k 2 



+ 



pik 2 
C2 
k 2 p 2 



iX B 



(3.1) 



/ 



We shall be interested primarily in the diffusional behavior contained within this matrix. This requires only the 
consideration of the limits given above. 



A. Expression for the structure factor 

As a consequence we can write the polymer-polymer structure factor for the case kR g <C 1. By recalling the 
definition of the matrix T, this must be the 1,1 element of the inverse, T _1 , 

G n (k, u) = 2 k 2 Pl (Cifc 4 + 2C 1 k 2 \ B p 2 + CiC 2 2 ^ 2 + CiAs V + Pi^fCiPi) L 2 

I (2 k 6 X B f 2 Pl L + k s + 2Lk 8 v Pl + A B 2 / 4 pi 2 fc 4 i 2 + Ci 2 w 2 fc 4 L 2 + Ci VC 2 2 ^ 2 
+X B 2 f 4 Pl 2 C 2 2 cu 2 L 2 + 2vk 4 Pl 2 X B 2 p 2 f 2 L 2 + 2vk 6 Pl 2 X B f 2 L 2 
+Qi 2 uj 2 X B 2 p 2 2 L 2 + v 2 k 4 Pl 2 ( 2 2 Lu 2 L 2 + 2 d 2 uj 2 k 2 X B p 2 L 2 
+2 C,xUJ 2 X B 2 p 2 f Pl Q 2 L 2 + 2 vk 2 Pl 2 ( 2 2 uj 2 X B f 2 L 2 + 2 v 2 k 6 Pl 2 X B p 2 L 2 
+v 2 k 4 Pl 2 X B 2 p 2 2 L 2 + v 2 k 8 Pl 2 L 2 + 2 k 2 ( 2 2 uj 2 X B f Px L + 2 k 4 X B 2 p 2 f 2 Pl L 
+k 4 ( 2 2 u 2 + 2 k 6 X B p 2 + k 4 X B 2 p 2 2 + 2 Lk 4 C 2 2 uj 2 v Pl 

+2Lk 4 X B 2 p 2 2 vp 1 +4Lk 6 vp 1 X B p 2 ) . (3.2) 

Similar results hold for the cases of the other structure factors of the system. This is the structure factor of the 
polymers in the coupled systems of ions and counterions, where there is an excluded volume interaction between the 
polymer segments, and a dynamic effective interaction, and dynamic alteration to the distribution of the polyion*. 
When the limits p 2 — > 0, X B — > are taken, the expression becomes equivalent to that of Fredrickson and Helfand.Ea 
Although this factor looks rather complicated, the form of the equation can be simplified significantly by noting 
the definition of the Debye parameter, 

k 2 = X B p 2 (3.3) 

and by introducing the quantity 

|- 2 = Pl v (3.4) 

which differs from the conventional definition of the screening parameter by Edwards0 by a factor of 4. 

A plot of this structure factor is shown in Figure 1, where convenient values have been used for the parameters. 
The polyelectrolyte peak is clearly visible and moves to greater values of k with increasing frequency. The peak has 
already been discussed in the context of statistical mechanics by various authors.u 

We find that broadening together with heightening in both k and ui directions of the polyelectrolyte peak is a result 
of either a decrease of the polyion charge density or of an increase of Debye parameter. 
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B. Diffusional behaviour 



By inspection of the fraction, equation 3.2, it is easy to identify that the numerator is a second order polynomial 
in u> and the denominator, when expanded is a polynomial of fourth order in the frequency. This suggests that the 
diffusion in the present approximation has two modes, which can be extracted by breaking the expression 3.2 into its 
partial fractions. The result of this manipulation is: 

G " = f4H + f4S> ( 3 - 5 ) 

1 i + 0T I 2 +U 

where the modes associated with the amplitudes Ai and A 2 have the following values which have been expanded to 
their lowest orders in k 2 : 

- - ,„2 _ PlP + P1P2V + L~ 1 p 2 
ClP2 + Ppi(2 



ryfc 2 ~ ' ' ; ^ = d c (3.6) 



r 2 ~-^f (/VC2 + C1P2). (3.7) 

This behavior to is to be anticipated in a system where there arc two components coupled by an interaction such as 
the Coulomb interaction, and has been previously deduced within the context of a far simpler assumption about the 
dynamics of the polyelectrolyte system. The mode T2 has already previously been identified as the plasmon mode. It 
is not proportional to k 2 at large length-scales and agrees with that found by other methods!] 

The cooperative diffusion_f;oefficient is given by the expression 3.6. Here we see that this is also given as in the 
paper of Vilgis and Borsalijfl with the limiting proportionality vp\ + 1/L. This is the expected behavior. We note 
that in the present scheme of approximation the cooperative diffusion coefficient, D c , depends upon the interaction 
only through the linear charge density of the chains and the countcrion concentration, whereas the plasmon mode is 
directly proportional to the Bjerrum length. 



IV. DISCUSSION 



In the first scheme of integration of section II, where the immediate integration over the fluid velocity field degrees 
of freedom leads directly to the nonlinear term in the polymer density, which is the only term coupling the nature of 
the interaction into the hydrodynamic description and vice-versa. Under the RPA for the transversal components too 
it is clear that this coupling term plays no contribution whatsoever, and then the dynamics of this system would be 
exactly as in the case derived by Fredrickson and HelfandEj For frequencies ujL 2 /tt 2 (i 3> 1 and at large length scales 
kR g < 1 a macroscopic solution shear viscosity is found 7/ = 77(1 + pi [it\r), with the intrinsic viscosity of the Rouse 
model [r)]n — L(il2rjE5 At length scales below the radius of gyration, Fredrickson and Helfand find a hydrodynamic 
screening length £ 2 oc ry/piCi? which screens the hydrodynamics such that at the smallest distances the viscosity is 
again 77. This is completely independent of the nature of the interaction between the polymer segments, as it contains 
only terms due to the connectivity of the polyion chains. 

However, the dynamic scattering intensity derived above can be implemented in an alternative computation of 
the reduced viscosity in the Rousej-as well as the Zimm cases by making use of the resulte-of Hess' and Klein's 
mode-mode coupling computations,^ as has been done by Vilgis, Borsali, and Benmouna,EIlE3 by implementation 

of the formula rj r = fo° dt dk <S l2 (k, t) grnm ' ^ e function H(k) represents the particle-particle 

correlation function. Hitherto this formula could only be implemented by estimating the first cumulant frequency 
for the dynamical scattering intensity, but here we have already computed this from the more complete dynamical 
formalism. We have also shown that in the appropriate limits, the behavior of our more complicated system is like 
that derived by more conventional means, so we expect the occurrence of a viscosity peak to be qualitatively retained 
in this formalism. 

It should be pointed out that the assumption of homogeneity of the separate constituents of charge and species in 
the solution, that these are in no way correlated when the Gaussian approximation of the collective coordinates is 
made, provides the approximation which leads to the Debye-Hiickel potential. This can also be derived in the other 
approximations, such as the linearization of the Poisson-Boltzmann equation. The RPA is valid for systems in which 
the temperature is high such that the various counterions can be understood to move almost independently and, and 
when the polymer density is high, such that the chains have a Gaussian end-to-end distribution. The success of this 
calculation has been tested in many and varied fields. As has been shown the RPA can be applied to the dynamical 
case too, with the complicating factor of the occurrence of additional fields; however, these are also quadratic. 
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In polyelectrolytes the effect of the condensation of counterions on the chains has been known for some time (see refs 
31, 32 and references therein). The critical onset of this condensation for monovalent ions which is derived from models 
of charged cylindrical objects occurs in Manning's theoryciUvhen Xgf > 1 (and at a lower clustering temperature in 
the case of solutions of rod-like polyelectrolytes; see LevinE3 and Levin and Barbosa.Eil) such that we consider here 
the case of high temperature and low charge density. Indeed, this other regime of counterion condensation on flexible 
objects is less well understood and would necessitate a different treatment than the RPA. 

Firstly, it is noticed that the presence of the hydrodynamic interactions generates a term to third order in the 
density of the fields. Attempts to approximate this as a linear or quadratic contribution to the integrand are bound to 
fail. This is because of the decoupling in the RPA of the fields relevant to the transverse and parallel k directions. The 
cubic term is responsible for the coupling of the hydrodynamic properties of the system to the remaining dynamic 
properties, and we can see how its neglect, makes the diffusive motion not feature any properties of the solvent, 
except through the definition of the friction parameter, £ = GTrrja, where a represents the particle radius. In order to 
understand more fully the role of the coupling (not only that of the connectivity) but especially of the nature of the 
interaction into the hydrodynamic properties, it becomes necessary to use an altered approach to the dynamics. This 
falls beyond the scope of this paper, and will be discussed in a subsequent publicationHj 

We have given the dynamical structure factor of the system under discussion, and shown that it reproduces the 
polyelectrolyte peak at sufficiently low frequencies. We have also discussed the diffusive behavior of the system, and 
reconfirmed the occurrence of the "plasmon" mode using the dynamical formalism. 
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APPENDIX: RPA RESULTS 

The RPA has been computed before by Fredrickson and Hclfand.El Here we present the relevant results for the case 

kR g < 1. 

2fc 2 Ci 

(pi(-k,-^)px(k,cj)) RPA = ^2 + fc4L -2 C 4 - 1 ) 

(p 2 (-k,-o;)p2(k,cj)) RPA = Jg + k ± ( 4 - 2 ) 

2k 2 f(i\ 

(p 1 (-k,- M )k. gl (k, M )) HPA %Ci+ , ta£ _ 1 (-J (4.3) 

2fc 2 U 2 
ujC,2 +ik 2 \ i 

Apart from the nontrivial remaining averages (7r a 7r a ) and ^7r a (l — kk)cr Q ) all remaining pairwise correlations are 
zero. 



( P2 (-k, - w )k • <x 2 (k, ^)) RPA = - ( H ) . (4.4) 
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FIG. 1. The dynamical structure factor for £1 = f 2 = 1, Pi = 1, / = 0.03, L — 100, v = 0, k = 0.01, and \ B = 1. The graphs 
show clearly how the peak in the dynamic structure factor broadens in the ^-direction. The maximum width in lj occurs at a 
momentum transfer vector larger than the location of the scattering peak. 
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